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Abstract 

The astrophysically important •^He(Q;,7)^Be reaction was studied at high 
^ \ energies where the available experimental data are in contradiction. A thin 

window ^He gas cell was used and the cross section was measured with the ac- 
^ ■ tivation method. The obtained cross sections at energies between Ec.m. = 1-5 

pH . and 2.5 MeV are compared with the available data and theoretical calcu- 

lations. The present results support the validity of the high energy cross 

section energy dependence observed by recent experiments. 

> 

O ■ Keywords: nuclear astrophysics, pp-chain, cross section measurement, 

j^ . '^He(Q;,7)^Be reaction 

O ' 1- Introduction 

m 

The '^He(a, 7)^Be reaction is an excellent example of those reactions which 
preserve their importance over several decades despite many experimental 
K^ \ and theoretical studies devoted to them. The crucial role played by the 

H ■ ^He(a,7)^Be reaction in solar hydrogen burning was realized more than 50 

ears ago, when H.D. Holmgren and R.L. Johnston measured its cross section 



yei 



l| and found a value being about a factor of 100 higher than estimated 
theoretically [2]. This result indicated that the 2nd and 3rd pp-chains |3| 
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Figure 1: Experimental data for the ^He(a;,7)^Be S-factor available at the end of the 20*'' 
century. The data are from [J, |5|, la, [Zl, |8|, lll| . The Gamow- windows relevant for the solar 
hydrogen burning via the pp-chain and for the big-bang nucleosynthesis are also shown at 
the bottom of the figure. 



can compete favorably with the first chain at the temperature of the solar 
core. 

Several experimental studies followed this pioneering work in the next 
decades aiming at the determination of the ^He(a, 7)^Be cross section at 
energies as low as possible in order to approach the astrophysically relevant 

This energy range, the so called 



V as possible in order 
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energy range 

Gamow window, lies between about 15 and 30keV at 15 MK temperature 
relevant for the solar core. The cross section at these energies is too small to 
be measured directly, theory-based extrapolations are therefore necessary to 
obtain the cross section at solar energies and hence the reaction rate. 

Figured] shows the summary of the experimental results available at the 
time when the first review paper on solar fusion reactions by E.G. Adel 
berger et al 
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was published. In order to avoid the strong energy de- 
pendence of the cross section caused by the Coulomb barrier penetration, 
the results are shown in the form of the astrophysical S-factoio. Based on 
these data, a recommended value of the zero energy extrapolated S-factor of 
S(0) = 0.53 ±0.05 keVb is quoted in Ref. 12|. The relative high uncertainty 



^For the definition of the S-factor see e.g. Ref. [3| page 172. 



of this value is on one hand due to the typically low precision of the individ- 
ual measurement at low energies and the large scatter of the data points (see 
Fig.[T]). On the other hand, the uncertainty is also increased by the fact that 
measurements carried out with two different experimental techniques (activa- 
tion and in-beam methods, see below) resulted in different S(0) extrapolated 



values 12 



At around the turn of the century the detection of solar neutrinos entered 
a precision era. The solar neutrino problem has been solved by the discovery 



of neutrino oscillation 13|, |lj] and the flux of ^Be and *B solar neutrinos 
can now be measured with a precision of about 3-4% by various neutrino 



detectors [15|, [ly, ll7|, ll8||. The measured neutrino fluxes can therefore be 
used to probe the solar core and test solar models if the underlying nuclear 
physics input, i.e. the cross section of pp-chain reactions is known with 
sufficient precision 19|]. Since the ^He(a, 7)^Be reaction is the starting point 
of the 2nd and 3rd chains from where the ^Be and ^B solar neutrinos originate, 
the uncertainty of its cross section must be reduced. 

The need for a more precise ^He(Q;,7)^Be cross section arose also from 
another field of nuclear astrophysics, namely the big-bang nucleosynthesis. 
Owing to the high precision measurement of some cosmological parameters, 
like the baryon-to-photon ratio of the universe carried out by the WMAP 
mission 
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primordial abundances of some light nuclei can be compared 
to observation with high accuracy. While good agreement is found between 
model predictions and observations for most of the nuclei produced in the 
big-bang, ^Li is significantly overproduced in the models. This so called ^Li 



problem is one of the biggest unresolved problems of big-bang theory [21 . 
At the observed values of the baryon-to-photon ratio of the universe the 
'^He(a,7)^Be reaction is responsible mainly for the production of ^Li. The 
^Be decays by electron capture to ^Li after the primordial nucleosynthesis 
ceasecl. The precise knowledge of the ^He(Q;,7)^Be reaction rate is crucial 
in order to exclude a solution of the ^Li problem based on this reaction rate. 
Answering the urgent call from these two fields of nuclear astrophysics, 
several new experiments studied the ^He(a,7)^Be reaction in the last 10 



years [22, l23|, M, M, |26|, |27|, l28| . The results of these "modern" ^B.e{a, jYBe 



experiments are shown in Fig.[2j Not surprisingly these experiments tried to 



^At lower values of the baryon-to-photon ratio ^Li would be mainly produced directly 
by the ■^H(Q;,7)^Li reaction. 




Nara Singh (2004) 
LUNA (2007) 
Brown (2007) 
ERNA (2009) 
Parker (1963) 
Kajino (1984) 
NoUctt (2001) 
Ncff (2011) 



ri; 



ITf*' 



-m 



1000 
Ec.m. 



/keV 



Figure 2: "Modern" experimental data for '^}ie{a,jyBe with total uncertainties. The 
data are from 22|, |23|, |2J, |25[_ 26l . 1271 . [28 1 . The results of three microscopic calculations 
by T. Kajino and A. ArimalsojTK.M. Nollett ^ and T. Neff ^ are also shown. The 
discrepancy between different experimental datasets and different calculations is clearly 
seen and gives the motivation of the present work. 



concentrate on the low energy range in order to be as close to astropliysical 
energies as possible. Only the ERNA collaboration using a novel approach 
of a recoil separator measured the cross section in a wider energy range up 
to Ec.m. = 3.1MeV [28|. Above about 1.3 MeV the only dataset previously 
available was that of P. Parker and R. Kavanagh j^. Fig. |2] also includes 
these data. As one can see, there is a clear disagreement between the two 
high energy datasets. The ERNA data are significantly higher than the 
P. Parker and R. Kavanagh data and the energy dependence of the S-factor 
is also different. 

This contradiction has a strong astrophysical consequence since the cross 
section must be determined at solar energies, well below the lowest available 
measurement and the extrapolations are based on theoretical calculations. 
The goodness of these theoretical models and thus the reliability of the ex- 
trapolation is assessed based on the comparison with experimental data at 
high energies where data are available. Moreover, the calculated S-factors are 
often fitted to available experimental data by a single rescaling parameter to 
obtain the extrapolated S-factor, even though such a rescaling does not have 
a strong physical justification for microscopic models 29|. As an example. 
Fig. [2] shows three theoretical cross section curves which are used to obtain 



cross sections at the Gamow window. The result of the microscopic models of 
T. Kajino and A. Arima (30|] and K.M. NoUett 3l(| were published when only 
the "old" experimental data were available. In lack of other data, their aim 
was to describe the data of P. Parker and R. Kavanagh at high energies. Af- 
ter the ERNA data became available, T. Neff carried out microscopic model 
calculations and was able to describe relatively well the "modern" experi- 
mental data 
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The zero energy extrapolated S-factors from the models 
differ by as much as 50 % (without any rescaling to experimental data). Sev- 
eral other theoretical models are available for the '^He(a, 7)^Be reaction (see 



e.g. Refs. |33|, M, M, M, M, l38|, |39|, |40|, lijj, |42|, |43|). The application of 



different theoretical models introduces a non-negligible uncertainty in the 
extrapolated S-factor. The new review paper on solar fusion reactions by 
E.G. Adelberger et al. J29| assigns a theoretical uncertainty of 4% to this 
value. In order to reduce this uncertainty, the high energy behavior of the 
'^He(Q;,7)^Be cross section must also be understood. 



2. The aim of the present work and the chosen method 

Presently there is no explanation to the difference between the two datasets 
at high energies. It is therefore necessary to carry out new experiments in the 
problematic energy range in order to provide more data for a better under- 
standing of the cross section behavior. The aim of the present work is thus to 
measure the cross section in the energy range between E^.m. = 1.5 MeV and 
2.5 MeV. 

There are two different ways to determine the '^He(a, 7)'^Be cross section 
which were both used several timeqj. The in-beam technique is based on the 
detection of prompt 7-radiation emitted when the fusion of the two Helium 
isotopes takes place in the direct capture process. The capture can lead 
to the ground or first excited states of ^Be and with the in-beam method 
partial cross section leading to these two states can be measured separately. 
The sum of the partial cross sections provides the total cross section. The 
other method is the activation which exploits the fact that the produced ^Be 
isotope is radioactive. It decays with a half-life of 53.22 ± 0.06 d [45] and the 



■^A third method is the direct counting of the created ^Be isotopes in a recoil sepa- 
rator. This technique was used up to now solely by the ERNA collaboration. A similar 
experiment is in preparation at the DRAGON facility at TRIUMF, Vancouver, Canada 



electron capture decay is followed by the emission of a single 7-radiation of 



478 keV with 10.44 ±0.04% relative intensity [45l|. The off-line detection of 
this 7-radiation can also be used to determine the cross section. This method 
provides directly the total cross section but gives no information about the 
branching ratio between the population of the ground state and first excited 
state. 

The set of "old" experiments showed an apparent discrepancy of about 
15 % between the results obtained with the two methods. The new exper 



iments of the LUNA collaboration 23, 24, 25, 27 and Brown et al. 26 



provided data with both methods and perfect agreement between the results 
were found. Therefore, in the present work only one technique was used. The 
activation method has some clear advantages over the in-beam technique, 
especially at higher energies. There is no need to determine the angular dis- 
tribution of the emitted 7-rays - as opposed to the in-beam method - and the 
beam induced background is of less importance. Therefore, the activation 
method has been chosen for the present work. 

Since the ^He(Q;,7)^Be reaction takes place between two isotopes of a 
gaseous element, the application of a gas target is necessary. At high energies, 
a thin window gas cell is a possible choice because the beam energy loss in 
the entrance foil is typically low compared to the beam energy and does not 
introduce a high uncertainty in the determination of the effective energy. A 
thin window gas cell is technically easier to handle and the determination of 
the number of target atoms can be done with high precision. The present 
study was thus carried out using a thin window gas cell. Further details 



about possible gas target designs can be found in [46 . 

3. Experimental details 

The experimental procedure used to determine the ^He(a,7)^Be cross 
section has been described in detail in a separate publication [46]. Here only 
the most important aspects of the experiments are summarized. 



The a-beam was provided by the MGC-20 cyclotron of ATOMKI [47 
After many tests, five irradiations were carried out between a-energies of 
4.0 and 6.3 MeV. The energy of the beam at the used beamline is defined 
to a precision of 0.3 %. The typical ^He"'""'" beam intensity was between 100 
and 250 pnA. The length of irradiations varied between 14 and 25 hours 
(see tableil]). Although the irradiation time was short compared to the half- 
life of the reaction product, the current integrator counts were recorded in 
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Table 1: Some parameters of the irradiations. 



Run no. 


Eq, 


Length of 


Total number 


Average target 




MeV 


irrad. [h 


of projectiles 


thickness [atoms/cm^ 


#1 


4.00 


24.3 


7.96-10^6 


2.16-10^9 


#2 


4.70 


24.7 


9.68-10^6 


2.33-10^9 


#3 


5.40 


13.1 


6.57-1016 


2.17-1019 


#4 


5.90 


16.2 


6.80-10i'5 


2.34-1019 


#5 


6.32 


22.7 


5.20-101'' 


2.35-1019 



multichannel scaling mode with 1 min time base in order to follow the changes 
in beam intensity. The current integrator was calibrated using a current 
generator and a precision current meter. 

The beam entered the target chamber through a beam defining aperture 
of 4 mm in diameter. After this aperture the whole chamber, including the 
gas cell, served as a Faraday cup in order to determine the number of pro- 
jectiles impinging on the target by charge measurement. A larger diameter 
aperture was placed behind the beam defining aperture and biased to -300 V 
to suppress the secondary electrons emitted from the target and aperture. A 
liquid nitrogen cooled trap was placed upstream the chamber to reduce the 
carbon build-up on the gas cell entrance foil. 

The ^He target gas was confined in a closed gas cell of 3 cm length. The 
^He isotopic enrichment was 99.95 %. The gas was filled into the cell before 
the irradiation and kept closed afterwards. The pressure of the gas was 
continuously measured with a type 722B MKS Baratron. The gas pressure 
was typically 300 mbar. The temperature of the gas cell was measured with a 
thermometer attached to the body of the cell. The back of the gas cell which 
served as the beam stop was directly water cooled. The typical temperature 
of the cell was about 25 °C determined by the cooling water temperature. 

The local density reduction of the gas when the beam is passing through, 
the so called beam heating effect, can alter significantly the effective number 
of target atoms. In the case of the present work the relatively low beam 
intensity and the low stopping power at the studied energies makes the beam 
heating effect rather low. The dissipated power of the beam was typically 
1 mW/mm. M. Marta et al. studied extensively the beam heating effect with 



a-beam penetrating a He target |48|. The obtained value for the necessary 



correction was found to be (0.0091 ±0.0019) (mW/mm) i. This translates 
into a 1 % beam heating correction in our case. Previously, J. Gorres et al. 



found a value of about a factor of two lower [49|. In the analysis of the 
present experimental data a beam heating correction factor of (1 ± 1) % was 
therefore assumed for all five runs. 

The beam enters the gas cell through a 1 /xm thick pinhole-free Ni entrance 
foil glued onto a stainless steel frame having an 8 mm diameter hole. The 
thickness of the foil was measured with ct-energy loss and perfect agreement 
with factory value of 1/im was found [46^. Tabled] lists some parameters of 
the irradiations. 

According to the reaction kinematics, the created ^Be nuclei move forward 
within a cone of a few degree opening angle. In the studied energy range and 
with the gas target thickness used, the ^Be nuclei have enough energy to reach 
the end of the cell and be implanted into the beam-stop. A removable Oxygen 
Free High Conductivity copper disk was used at the beam stop as catcher 
to collect the created reaction products. The catcher material was chosen 
based on several test measurements aiming at finding parasitic production 
of ^Be in them. For details, see [46|. GEANT4 simulations proved that 
more than 99.7% of the produced ^Be reaches the catcher within a circular 
spot of 24 mm in diameter around the beam axis. At the studied energies 
backscattering loss of ^Be from the catcher is well below 1 % and it is therefore 



neglected [46 



After the irradiations the Cu catcher was removed from the chamber and 
the collected ^Be activity was measured with a 100 % relative efficiency HPGe 
detector equipped with a complete An low background shielding. In order to 
maximize the detector efficiency, the catcher were placed in a close geometry 
in front of the detector. The distance between the catcher and the detector 
end cap was 1 cm. Since ^Be is a single line 7-emitter, no true coincidence 
summing effect is present. In order to avoid coincidence summing also in 
the case of calibration sources, the absolute efficiency of the 7-detector was 
measured with single line sources only. These sources were the following: 
^"^Mn, ^^'^Cs, ^^Zn and ^Be. Fig.[3] shows the measured efficiencies with these 
sources and the fitted curve. In order to take into account the different 
spatial distribution of the calibration sources and the actual ^Be source on 
the catcher, GEANT4 simulations were carried out. The final efficiency at 
478 keV was determined to be (8.87 ±0.44) %. 

Table[2]lists some parameters of the 7-countings. Fig.lHshows the relevant 
part of the 7-spectra of the five runs as well as the laboratory background 
spectrum. The spectra are normalized to counting time. 



t6 0.07 







E^ I keV 

Figure 3: Absolute efE.ciency of the HPGe detector as a function of energy. In order to 
avoid true coincidence summing, only single line calibration sources were used. The good 
log-log straight line fit to the measured data strengthens the reliability of the measured 
efficiency at 478 keV with the calibrated ^Be source. 



Table 2: Some parameters of the 7-countings. The waiting time is the time elapsed 
between the end of the irradiation and the start of the counting. 



Run no. 


Waiting 


Counting 


^Be net 




time [h 


time [h 


area [cts 


#1 


1 


192 


4559 ±178 


#2 


106 


168 


6422 ±216 


#3 


8 


288 


8630 ±236 


#4 


81 


192 


6681 ±239 


#5 


63 


168 


5661 ±175 
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Figure 4: Relevant part of the measured 7-spectra in the five runs. In the lower right 
panel the laboratory background is plotted. The spectra are normalized to counting time. 
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The cross section a was calculated using the following formula: 



where T [atoms/cm^] is the target thickness, [1/s] is the beam intensity, 
A is the ^Be decay constant, e^ is the gamma detection efficiency at 478 keV 
and r] is the decay branching to the ffist excited state in ^Li leading to the 
478 keV 7-radiation. A is the measured 478 keV net peak area. The time 
constants ti, ty^, t^ are the length of the irradiation, the waiting time between 
irradiation and counting and the length of the counting, respectively. 

4. Results and discussion 

Table[3] summarizes the experimental results. The second column shows 
the effective center-of-mass energies. The effective energies were calculated 
based on the following steps. The energy loss in the Ni entrance foil obtained 
from the SRIM code [50|] was subtracted from the beam energy. The thickness 
of the foil was taken to be 1.01 ±0.04//m [46]. Then the energy loss in the 
target gas was considered. Since the energy loss in the gas is rather low 
(typically around 150 keV) and the S-factor does not change much within this 
energy interval, a constant S-factor over the target thickness was assumed 
and the effective energy was calculated based on the procedure given e.g 
in chapter 4.8 of |3|. Finally, the effective energy was converted into the 
center-of-mass frame. The uncertainty of the effective energy is obtained by 
assigning 0.3 % uncertainty to the primary beam energy, taking into account 
the energy straggling and using 3.9% uncertainty for the stopping power 



values [Ml- In principle, the energy straggling cannot be treated the same 
way as the other uncertainties since it represents only the broadening of the 
beam energy distribution. In order to quote a conservative estimate of the 
uncertainty of energy, and taking into account its small contribution to the 
final uncertainty, the value of the energy straggling is, however, included in 
the calculation of the energy uncertainty. 

The uncertainty of the cross section values is calculated as the quadratic 
sum of the statistical error of the ^Be peak area and several systematic un- 
certainties. These uncertainties are listed in TableHl In Table[3]the statisti- 
cal and total systematic uncertainties, respectively, are quoted separately in 
parenthesis. The uncertainty of the S-factor values also includes the effective 
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Table 3: Results of the present work in the form of cross section and S-factor. The statis- 
tical and total systematic uncertainties, respectively, are quoted separately in parentheses. 
See text for fu rther details. 



Run no. 



Cross section 
[keV] [/ibarn] 



S-factor 
[eV barn] 



#1 


1473 ±15 


2.96(12)(17) 


313(12)(19) 


#2 


1791 ±14 


3.78(13)(22) 


327(11)(19) 


#3 


2115 ±14 


4.68(13)(28) 


351(10)(21) 


#4 


2338 ± 14 


5.08(18)(30) 


354(13)(21) 


#5 


2527 ±14 


6.06(19)(36) 


401(12)(24) 



Table 4: List of uncertainties on the cross section (S-fa 


ctor) determination 


Type Source of 
uncertainty 


Value in 
percentage 


statistical ^Be peak area 


2.7-3.9% 



7-detection efficiency incl. source distribution 4 % 

number of target atoms* 3 % 

systematic current integration 3 % 

beam heating 1 % 

effective energy** 0.4-1.1 ' 



total systematic 



5.9 



* includes gas pressure, temperature and gas cell length uncertainties 
** affecting only the S-factor, not the cross section 



energy uncertainty, however, since the effective energy is rather well defined, 
this source of uncertainty is small compared with the others listed in tablelH 

FigureO shows the results in the form of S-factor together with the data 
of the "modern" experiments and that of P. Parker and R. Kavanagh. Along 
with the experimental data, the calculated S-factor of the three microscopic 
models introduced in Fig. [2] are also plotted without any rescaling. 

When the present study was in progress, M. Carmona-Gallardo et al. 
published their results on '^He(a,7)^Be in a similar energy range providing 



cross section values at three energies [52[. The motivation of that work was 



similar to ours. The results of this new experiment are also included in Fig.|5l 

The present data support the findings of the two recent high energy 

datasets by the ERNA collaboration [28| and M. Carmona-Gallardo et al. 



52| | . The rather low values and fiat energy dependence of the P. Parker and 
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Figure 5: Results of the present work compared with "modern" data and some theoret- 
ical calculations. Only the statistical uncertainties are plotted with the exception of the 
P. Parker and R. Kavanagh data where different sources of uncertainty are not quoted 
separately. See text for discussion. 

R. Kavanagh data are not reproduced by the new experiments. 

In summary, the results presented in this paper give a significant con- 
tribution to the better understanding of the high energy behavior of the 
'^He(a,7)^Be cross section. In the review paper of E.G. Adelberger, et al. 
|29[ the available ^He(Q;, 7)^Be data have been critically reviewed, but the S- 
factor fitting was restricted to energies below 1 MeV. A global analysis of the 
experimental data also at higher energies and their comparison with theory 
would also be needed for a more complete description. Such an analysis is 
beyond the scope of the present work. 
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